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Thin film magnetic heterostructures with competing interfacial coupling and Zeeman energy pro¬ 
vide a fertile ground to study phase transition between different equilibrium states as a function 
of external magnetic held and temperature. A rare-earth (RE) / transition metal (TM) ferro¬ 
magnetic multilayer is a classic example where the magnetic state is determined by a competition 
between the Zeeman energy and antiferromagnetic interfacial exchange coupling energy. Techno¬ 
logically, such structures offer the possibility to engineer the macroscopic magnetic response by 
tuning the microscopic interactions between the layers. We have performed an exhaustive study 
of a nickel/gadolinium system by using the element-specihc measurement technique x-ray magnetic 
circular dichroism, and determined the full magnetic state diagrams as a function of temperature 
and magnetic layer thickness. We explain our result based on a modihed Stoner-Wohlfarth formal¬ 
ism and provide evidence of a thickness-dependent phase transition to a magnetic fan state which 
is critical in understanding magnetoresistance effects in RE/TM systems. The results provide im¬ 
portant insight for spintronics and superconducting spintronics where engineering tunable magnetic 
inhomogeneity is key for certain applications. 

Modified magnetic interaction at the interfaces be¬ 
tween different ferromagnetic materials can be utilised 
to engineer materials with magnetic properties that are 
significantly different from the bulk. A classical exam¬ 
ple of such a system is a rare earth / transition-metal 
(RE/TM) multilayer since in the presence of an external 
magnetic field the layers favor parallel alignment but the 
antiferromagnetic coupling between the layers can lead 
to novel magnetic states and phase transitions between 
them. Ferrimagnetic alloys and multilayers (favored an¬ 
tiparallel alignment between layers) are of great inter¬ 
est from both a fundamental and an applied point of 
view. Ferrimagnetic alloys were until very recently the 
only magnetic system that showed ultrafast magnetiza¬ 
tion switching induced by a femtosecond laser pulse [1-3], 
and engineered magnetic multilayers are key components 
of superconducting spin valves [4-6]. 

The theoretical groundwork of RE/TM systems was 
established by Camley and Tilley [7, 8]; they modelled 
RE/TM multilayers and found that due to the differ¬ 
ent Curie temperatures (T/) of the two magnetic com¬ 
ponents several phases occur as the effects of tempera¬ 
ture and applied magnetic field are incorporated. At low 
temperatures (T < Tc{RE)), the RE aligns with the ap¬ 
plied field and the TM aligns antiparallel to the RE and 
the field (RE-aligned state), while at higher temperatures 
(T > Tc{RE)) the magnetization of the RE decreases 
and a second-order phase transition to the TM-aligned 
state occurs. This work was later confirmed experimen¬ 
tally [9, 10], however, Camley also predicted that at the 


interface between bulk Fe and a thin layer of Gd (five 
atomic layers), a twisted magnetic state could form at 
some values of applied field [11], although this is still an 
open question. 

Due to the interface-driven nature of RE/TM systems 
it is difficult to measure the exact magnetic structure 
of a system with bulk magnetic measurement techniques 
such as vibrating sample magnetometry and magnetore¬ 
sistance measurements [12-14]. Additionally, at low tem¬ 
peratures the magnetization of a RE layer in a typi¬ 
cal RE/TM bilayer will be far larger than that of the 
TM layer, making it extremely difficult to deconvolute 
the two signals. Therefore, an element specific measure¬ 
ment technique such as x-ray magnetic circular dichroism 
(XMCD) is ideally suited to measuring such a system. 
The energy of the x-rays used can be tuned to the absorp¬ 
tion edges of particular elements, meaning only the mag¬ 
netic response of specific layers is measured at any one 
time. Previous studies using XMCD have mainly con¬ 
firmed the predictions of Camley and Tilley [7, 8]: Barth 
et al. observed evidence of the transition between the 
RE- and TM-aligned state as a function of temperature 
in a Ni/Gd bilayer (7.5 nm and 5 nm thick, respectively) 
[10], while Koizumi et al. observed the same behavior in 
a Gd/Fe (2 nm/2 nm) 5 o superlattice measured at 20 K 
and room temperature [15]. 

In the present study we present extensive XMCD re¬ 
sults of Ni/Gd/Ni films allowing us to fully map the mag¬ 
netic states as a function of applied field and temperature 
for a wide range of different Ni and Gd thicknesses (Figs 
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1(a) and (b)). Our results demonstrate the full richness of 
the magnetic phase diagram of RE/TM systems and shed 
new light on previous work that interpreted results based 
on an incomplete picture of RE/TM behavior [12, 16]. 

We measured 15 samples by XMCD, each with a differ¬ 
ent layer thicknesses, throughout a range of temperatures 
between 6 K and room temperature [17]. The measured 
samples displayed one of three behaviors as illustrated 
in Fig. 1(c) which shows element-specific magnetiza¬ 
tion versus applied field {M{H)) hysteresis loops at high 
and low temperatures for three different samples that are 
characteristic of each behavior. In the first behavior the 
Ni layers follow the applied magnetic field (positive mag¬ 
netization in positive field) while the Gd layer remains 
antiparallel to the Ni and the magnetic field. This does 
not change as the temperature is increased. The hystere¬ 
sis loops shown are from a sample with total Ni thickness 
of 5 nm, and Gd thickness of 4.5 nm [18]. In the second 
behavior, the opposite is observed at lower temperatures: 
the Ni is antiparallel to the field while the Gd is parallel. 
Above a transition temperature, this state is reversed to 
the same as that in behavior 1. The thicknesses of the 
Ni and Gd in this sample are 6 nm and 8 nm, respec¬ 
tively. These two behaviors can be understood in terms 
of the variation in the influence of the interfacial layers 
in the multilayer compared to that of the bulk. In the 
thinnest samples, the interfacial coupling will dominate 
and if the ratio between the Ni and Gd thicknesses favors 
the Ni, the Ni will control the behavior of the Gd despite 
the higher magnetic moment of the Gd at lower tem¬ 
peratures. For slightly thicker Gd layers, although the 
interfacial coupling still dominates the behavior of the 
system, the Gd will dominate the Ni while its magnetic 
moment is not suppressed by the increasing temperature. 

The third behavior is more complex, with the Ni M(H) 
loops at low temperature displaying a distinctive “N”- 
shape, indicative of significant magnetic inhomogeneity 
in Ni during magnetization reversal. Meanwhile, the Gd 
layer undergoes a gradual transition over a temperature 
range of tens of kelvins from aligned parallel to the Ni 
to antiparallel. The sample displaying this behavior in 
the figure had Ni and Gd thicknesses of 8 nm and 6 nm, 
respectively, and is the same sample as that for which 
the XMCD spectra are shown in Fig. 1(b). We do not 
believe this behavior in the Ni has been observed before; 
the magnetization crosses zero three times on each branch 
of the hysteresis loop. 

In order to understand the micromagnetic structure 
that corresponds to these hysteresis loops, we have mod¬ 
elled the TM/RE/TM multilayer. The model uses the 
Stoner-Wohlfarth framework by Mauri et al. [19], and 
that of Camley and Tilley, wherein an iterative process 
is used to find the ground state of a stack of Ni and Gd 
atomic layers. However, our model includes a term de¬ 
scribing the magneto-crystalline anisotropy energy since 
some previous work had not been able to use Camley 


and Tilley’s model to explain the behavior of an RE/TM 
multilayered system [20], citing the possibility that only 
including the Zeeman and exchange energy terms leads 
to limitations in explanatory power. We solve the follow¬ 
ing Hamiltonian to find the angle of one spin at which 
the energy is a minimum: 
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where g is the Lande g-factor, /xb is the Bohr magne¬ 
ton, ^0 is the permeability of free space, M is the average 
magnetization of the layer (in our case M is equal to the 
value of the magnetic moment of the atom for a perfectly 
ordered material), </> is the angle of a magnetic moment 
with respect to the applied magnetic field direction, K is 
the magnetocrystalline anisotropy energy, 0 is the direc¬ 
tion of the magnetic easy axis, n is the number density 
of the material of the layer, J is the exchange constant 
of the layer, and S is the magnitude of the magnetic mo¬ 
ment of an atom in units of the Bohr magneton. We then 
repeat the process for each spin in the Ni/Gd/Ni stack 
until successive iterations over the whole stack differ by 
no more than a factor of l/IOOO. A detailed description 
of the model can be found in the Supplemental Material. 
Figure 2 shows a comparison between the results of the 
model and the temperature dependence of the third type 
of behavior, which match very closely despite the simplic¬ 
ity of the model. Although the model can reproduce all 
three of the behaviors discussed above, here we focus only 
on the thrid type. It should be noted that one of the few 
assumptions made for the model is that of strict antifer¬ 
romagnetic coupling (negative exchange energy) between 
the Ni and Gd layers. The success of this assumption con¬ 
firms that no more complicated effects are influencing the 
behavior of the interface in order to understand the com¬ 
plex M{H) loops. We can see that at low temperatures, 
the Ni magnetization increases with increasing field, but 
abruptly reverses as the field crosses zero, before saturat¬ 
ing. This behavior is reproduced by the model, as is the 
thermal evolution of the hysteresis loops. As the tem¬ 
perature increases, the “N”-shape opens up until at ~ 
30 K the Ni loop looks more conventional, except with 
inflection points on the loop at the same field value as 
the inflection points at low temperatures. 

In the Gd measurements, we can now see the details of 
the transition from aligning parallel to aligning antipar¬ 
allel with the field. The saturation magnetization. Mg, 
at negative field increases while the opposite is true for 
the Mg at positive fields. Interestingly, for most of the 
temperature range shown the saturation magnetization is 
not the highest value of magnetization within the loop; at 
saturation there must be significant Gd moments canted 
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away from the applied field direction. This feature, com¬ 
bined with the unusual shape of the Ni loops leads us 
to conclude that some form of magnetic inhomogeneity 
exists in the Ni/Gd/Ni structure. 

The model we use allows us to extract the micro- 
magnetic structure from the system that produced the 
hysteresis loops since the modelled system is a one¬ 
dimensional stack of magnetic moments. We assume that 
the moments in each layer are parallel, and so can ap¬ 
ply a normalisation factor to each layer’s magnetization. 
This helps to reduce the computational complexity of 
the model, and since the measured samples are polycrys¬ 
talline, film orientation effects will be negligible. The 
model does not take into account the crystal structure 
since in a polycrystalline film orientation effects will be 
negligible. The sum of the projections of the magnetic 
moments of each of the elements onto the direction of the 
applied field gives the total magnetization at each field 
value, as shown in the hysteresis loops. But the original 
angles of the each magnetic moment can also be plot¬ 
ted individually to give vector diagrams that show the 
inhomogeneity in the system at particular field values. 

These vector diagrams are shown in Fig. 2(b). They 
show the angles of the magnetization for each atomic 
layer calculated by the model at a particular field value 
and temperature; the field values are indicated on the 
corresponding hysteresis loops. From the vector dia¬ 
grams we can see that the reason for the shape of the 
Ni hysteresis loops at low temperature is the effect the 
Ni/Gd interface has on the Ni layers. At saturation fields 
the Ni magnetization is not at its maximum value be¬ 
cause the moments at the interface have been canted 
away from the field direction due to the competition 
between the inter-, intra-layer exchange energiesand the 
Zeeman energy. The sharp decrease in magnetization in 
the Ni at small positive fields is due to the Gd exerting its 
influence on the Ni when the Gd reverses. The reversal 
of the dominant Gd “drags” the Ni round to remain an¬ 
tiparallel to the field before the increasing Zeeman energy 
of the Ni at increasing fields forces the degree of winding 
near the interface to be suppressed. By examining the 
model in this way, a unique understanding of the system 
and its magnetic microstructure can be gained. 

As well as yielding information about the magnetic mi- 
crostucture of the system, the model can be used to ex¬ 
tract other magnetic properties. For example, one fitting 
parameter in the model is the Tq of the Gd. Due to the 
large difference in the Curie temperatures between RE 
and TM ferromagnets, and the temperatures at which 
the experiment was performed, the value of the Ni sat¬ 
uration magnetization as a function of temperature in 
the model is presumed to be constant. However, since 
the Tc of bulk Gd is very close to room temperature, 
the decrease in magnetization as a function of tempera¬ 
ture must be taken into account in the model, and the 
effective Tc of Gd will be further suppressed due to the 


geometry of the Gd layer. The value of Tc can be ex¬ 
tracted from the model; Figure 3(a) shows the agreement 
between the experiment and the model when the Tc of 
Gd in the model is set to 35 K. The suppression of the 
Curie temperatures of ferromagnets as a function of film 
thickness is well documented [21-23]. The AMg values 
were calculated by averaging the maximum and mini¬ 
mum magnetization values within each loop and taking 
the difference between them. It should be noted that 
the magnetization values have been normalised to be¬ 
tween zero and one, and that therefore the rise in the 
Ni magnetization with temperature is not comparable to 
the decrease of Gd magnetization. In fact, the rise in the 
Ni originates due to the method of measurement; XMCD 
only detects the projection of a sample’s magnetization 
onto the k-vector of the x-ray, so as more Ni moments 
are able to align with the applied field ,which is parallel 
to the x-ray’s k-vector, the measured magnetization will 
increase. 

Figure 3(b) shows a magnetic state diagram for a par¬ 
ticular sample, that can be extracted from measured hys¬ 
teresis loops due to their element-specific nature. The Ni 
and Gd hysteresis loops are first normalized, so the max¬ 
imum of the modulus of the difference is 2 (AMg = 2); 
this corresponds to the Ni and Gd having oppositely sat¬ 
urated magnetizations, while AMg = 0 corresponds to 
parallel Ni and Gd. We know that when the Ni and Gd 
layers are antiparallel there is no inhomogeneity at the 
interfaces, but the inhomogeneity will nucleate at the in¬ 
terfaces as the angle between the layers changes. Mag¬ 
netic state diagrams such as these will be invaluable in 
the design of devices in a variety of fields that require 
engineered magnetic inhomogeneity. 

We must point out that in all the Figures showing the 
output of the model, the magnetic field values that can 
be compared with those from the experiment, coercive 
field values for instance, are larger than their counter¬ 
parts by an approximately constant factor of 4. This 
difference arises due to the fact that simple models do 
not include nucleation points present in physical samples 
that facilitate magnetization reversal at lower magnetic 
field values [24-27]. This discrepancy should not affect 
the results of the model in any way other than a scaling 
of the equivalent fields. 

The measured samples were grown in four different 
growth runs, but we are confident that the observed ef¬ 
fects are due to thickness dependence, and not stochastic 
variations in the growth conditions. Firstly, because sam¬ 
ples from different growth runs display the same types of 
behaviors. Secondly, imaging the magnetic landscape on 
a micron-scale for three samples using a scanning SQUID 
microscope [28], a sensitive probe of magnetic flux, shows 
consistent resolution-limited domain structure and uni¬ 
form magnetization strength for a variety of location in 
each sample (see Supplemental Material for more infor¬ 
mation). 
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In summary, we have performed a comprehensive ex¬ 
perimental and theoretical study of the RE/TM system, 
using both a more sophisticated model than has previ¬ 
ously been utilised and an extremely powerful measure¬ 
ment technique (XMCD). The use of the model allows 
a detailed view of the micromagnetic structure that is 
responsible for any of the three behaviors that were ob¬ 
served. The strict antiferromagnetic exchange coupling 
between RE and TM ferromagnets and its competition 
with the intralayer exchange energy and the Zeeman en¬ 
ergy has been shown to lead to significant inhomogene¬ 
ity in the system. The relative strengths of these ener¬ 
gies, and thus the behavior of the system can be tuned 
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FIG. 1. (a) Diagram of experimental setup showing the struc¬ 
ture of the samples, the direction of the incident x-rays, and 
the position of the measuring photodiodes (7o and I), (b) 
Examples of XMCD spectra taken at ~6 K for Ni (i) and Gd 
(ii). The XMGD was measured at the Ni L 3 and Gd M 5 ab¬ 
sorption edges, respectively, marked by dotted lines. The two 
signals and fj.~ are from positive and negative saturation, 
while Ajj, is the difference between the two, vertically offset 
for clarity, (c) Hysteresis loops measured between 5 and 10 K 
(low T), and 45 and 65 K (high T), that are representative of 
the three observed behaviors. The (total) Ni and Gd thick¬ 
nesses of each of the samples are (from left to right): 5 nm 
Ni and 4.5 nm Gd, 6 nm Ni and 8 nm Gd, and 8 nm Ni and 6 
nm Gd. The last sample is the same as that from which the 
XMGD spectra are shown in (b). 
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FIG. 2. Comparison between the hysteresis loops measured by XMCD and those produced by the model, as well as detailed 
view of the micromagnetic structure of the sample at certain temperatures and field values extracted from the model, (a, c) 
Hysteresis loops, experimental and theoretical for Ni and Gd, respectively. The magnetization has been normalised by the 
highest and lowest values of magnetization within the loop to ease comparison. Shaped markers on some theoretical loops show 
the temperature and field values of the vector diagrams in (b). (b) Vector diagrams showing the direction of magnetization of 
each atomic layer in the model system used to produce the hysteresis loops. The magnitude of the magnetization is not shown. 
Black arrows represent Ni layers, red Gd, as labelled in the top left. The experimental hysteresis loops are from the same 
sample as that shown in Fig. 1(b) and (c) (behavior 3), and has (total) Ni and Gd thicknesses of 8 nm and 6 nm, respectively. 
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FIG. 3. (a) Difference between the average maximum and 
minimum magnetizations in experimental and theoretical hys¬ 
teresis loops. The connected dots represent experimental re¬ 
sults, while solid lines represent the results of the model. The 
plot is generated by taking half of the difference between the 
average maximum and minimum values of magnetization in a 
hysteresis loop at a particular temperature. The whole plot 
has been normalised to its own maximum and minimum for 
ease of comparison, (b) Phase diagram showing the degree 
of inhomogeneity present in an experimental multilayer as a 
function of applied field and temperature. The phase diagram 
is constructed from experimental hysteresis loops by taking 
the difference between the normalised magnetization of Ni 
and Gd at each field value in a loop. This is then repeated 
for every temperature. 










